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Abstract 

We compare the theoretical status and the numerical predictions of two approaches 
for heavy quark production in the high energy hadron collisions, namely the conven- 
tional LO parton model with collinear approximation and fcr-factorization approach. 
The main assumptions used in the calculations are discussed. To extract the differ- 
ences coming from the matrix elements we use very simple gluon structure function 
and fixed coupling. It is shown that the fcT-factorization approach calculated formally 
in LO and with Sudakov form factor accounts for many contributions related usually 
to NLO (and even NNLO) processes of the conventional parton model 
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1 Introduction 



The investigation of the heavy quark production in high energy hadron colhsions provides 
a method for studying the internal structure of hadrons. Reahstic estimates of the cross 
section of the heavy quark production, as well as their correlations, are necessary in order 
to plan the experiments on existing and future accelerators. 

QCD is the theory of quark and gluon interactions. The description of hard hadron 
collisions is possible only at the phenomenological level, reducing the hadron interactions to 
the parton-parton ones, treated in terms of the hadron structure functions. The hard process 
cross section is expressed as the convolution of the parton distributions in the colliding 
hadrons with the cross section of the elementary sub-process given by the square of the 
matrix element, calculated in the perturbative QCD. Such an approach can be justified only 
in the Leading Logarithmic Approximation (LLA) in the leading order (LO), or in the next 
to leading order (NLO) and so on. Strictly specking, it is in disagreement with the quantum 
field theory, where the amplitudes, rather than the probabilities, should be added and/or 
multiplyed. Thus all theoretical approaches to the high energy hadron interactions are in 
some sense phenomenological ones. 

There are several phenomenological approaches to the calculations of hard processes cross 
sections in hadronic reactions. 

The most popular and technically simplest approach is the so-called QCD collinear ap- 
proximation, or parton model (PM). In this model all particles involved are assumed to be 
on mass shell, having only longitudinal components of the momenta, and the cross section is 
averaged over two transverse polarizations of the incident gluons. The virtualities of the 
initial partons are taken into account only through their densities (structure functions). The 
latter are calculated in LLA using DGLAP evolution equation and available experimental 
data. The probabilistic picture of "frozen" partons underlies this way of proceeding. The 
cross sections of QCD subprocess are calculated usually in the LO, as well as in the NLO 

1^, 1^, ^, 1^. However the transverse momenta of the incident partons are neglected in the 
QCD matrix elements in the direct analogy with the Weizsaecker- Williams approximation in 
QED. It allows to describe quite reasonably the experimental data on the total cross sections 
and one-particle distributions of produced heavy flavours, however it can not |^ repro- 
duce, say, the azimuthal correlations |§] of two heavy quarks, as well as the distributions 
over the total transverse momentum of heavy quarks pairs which are just determined by 
the transverse momenta of incident partons. 

The simplest way to incorporate the transverse momenta of the incident partons for the 
proper description of heavy quark correlations was suggested in 0, as the random shift of 
these momenta according to certain exponential distribution. However, this procedure has 
no serious theoretical background. While the shift of the order of Aqcd ((^t) ~ 300 MeV) 
looks to be reasonable as having an origin in confinement forces at large distances, the value 
(kx) ~ 1 GeV or more should be explained in terms of the perturbative QCD. Moreover, 
the averaged value of the shift seems to be dependent on the initial energy, kinematical 
region, etc. The simultaneous description of one-particle distributions over the transverse 
and longitudinal momenta seems to be impossible ||10|| . 
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Another natural method to account for the incident parton transverse momenta is referred 



to as fcy-factorization [T^, [T2|, |T3|, |14|, |T5|, or the theory of semihard interactions ||16| , |T8 
|19|. Here the Feynman diagrams are evaluated for the virtual incoming partons with all 
possible polarizations. Since there are no reason to neglect the gluon transverse momenta 
qiT and q2T in comparison with the quark masses and transverse momenta in the small x 
domain, the QCD matrix elements of the sub-processes become more complicated. We obtain 
them in the LO, the NLO calculations being much more technically hard in this approach 
than in the conventional PM. On the other hand, the essential and, most probably, the 
major part of the NLO (and part of NNLO) corrections to the LO PM related to the finite 
values of parton transverse momenta are already included into LO contribution in case of 
fcr-factorization. Thus one can hope that NLO, NNLO... corrections will not be too large 
in /c^-factorization, making it to be more constructive. 

The matrix element describing collision of the off-shell partons is not, of course, gauge 
invariant. This is the shortness of /cT-factorization. One can use it only in a "physical" 
(axial or planar) gauge, where the LO Feynman graphs responsible for the parton (DGLAP 
or BFKL) evolution take a ladder-like form, and the hard cross section can be factorized, that 
is written as a convolution of the initial parton wave functions with the hard sub-process 
matrix element squared. Recall, that an unintegrated parton distribution in this case is 
given not just by the derivative of the integrated parton density with respect to the scale, 
but includes the double logarithmic (Sudakov) form factor as well. This point has been 
discussed in more detail in pO[ (see also Sect. 3 of the present paper). In Feynman gauge 
this double logarithmic contribution comes from the graphs where the soft gluons embrace 
the "hard" blob, as it is shown, for example, in Fig. lb (contrary to the PM case shown in 
Fig. la). 

Note, that when the initial energy grows, the momentum fraction x carrying by the "ac- 
tive" (participating in hard interaction) parton decreases. As a rule the smaller x the larger 
number of steps the parton passes in the evolution. Having emitted a lot of evolution gluons 
the parton gains a large transverse momentum fc^. It explains why a large phenomenological 
kT shift is needed to describe the high energy data in terms of the PM. 

Today it is clear that the fcr-factorization approach is self-consistent and predicts the 
results in reasonable agreement both with the experimental data and with the collinear 
approximation when the last one can be applied. However, the detailed numerical analysis 
of hard processes in the framework of /cr-factorization is not provided yet. 

The predictions of all phenomenological approaches depend significantly on the quark 
and gluon structure functions. The last ones are more or less known experimentally from 
the data of HERA, but unknown at very small values of Bjorken variable x < 10~^. However 
it is just the region that dominates in the heavy quark production at high energies^. 

The main goal of this paper is to compare the results of the conventional parton model 
and the fcr-factorization approach and to explain the source for the differences in their 
predictions. We present the main formalism of these approaches as well as the numerical 



^For example, in the case of charm production, rric = 1.4GeV, at LHC, ^/s = 14 TeV, the product xiX2 
of two gluons (both xi and X2 are the integral variables) is equal to 4 • 10"^ and the applicapability of the 
existing structure functions at so small x is not clear, see discussion in 
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results. To clarify the results we consider very high energies, including the non-realistical 
one = 1000 TeV, and use the simple "toy" gluon structure function of the nucleon. 



2 Conventional parton model approach 



The conventional parton model expression for the heavy quark hadropro duct ion cross sec- 



tions has the factorized form 22 



da{ab^QQ) = j dxidxjGa/i{xi, HF)Gb/j{xj, iJ,F)d a{ij ^ QQ) , (1) 

where Ga/i{xi, fip) and Gb/j{xj, fip) are the structure functions of partons i and j in the 
colliding hadrons a and b, fip is the factorization scale (i.e. virtualities of incident partons) 
and da{ij — >■ QQ) is the cross section of the subprocess calculated in the perturbative QCD. 
It can be written as a sum of LO and NLO contributions, 



daizj^QQ) = :^i^if^f + A7ra,{f^nm> + n;>\n{^^^^^ (2) 

where /x/j is the renormalization scale, and the functions fj;°\ fj;^ as well as f-p actually 
depend only on a single variable 

ArriQ 

P ~ Z ; = XiXjSab ■ (3) 

The expression (|l]) corresponds to the process shown schematically in Fig. la. The main 
contribution to the heavy flavour production cross section at small x is known to come from 
gluons, i = j = g. 

The principal uncertainties of any numerical QCD calculation are the consequences of 
unknown scales 0, /Ui? and /x/j and the heavy quark mass, mg. Both scales (sometimes 
they are assumed to be equal) are to be of the order of hardness of the treated process, 
however which value is better to take, mg, ttit = + or s, remains to be unknown. 

Principally, the uncertainties should not be essential because of the logarithmical dependence 
on these parameters. Unfortunately the masses of c and b quarks are not large enough and it 
leads to numerically large uncertainties (see e.g. 0) at the existing energies for fixed target 
experiments. 

Another principal problem of parton model is the coUinear approximation. The transverse 
momenta of the incident partons, QiT and Qjt are assumed to be zero, as is shown in Fig. la, 
and their virtualities enter only via structure functions; the cross section daiij — > QQ) is 

^These uncertainties have to disappear after all high order contributions are summed up. There is the 
opinion that the strong (weak) scale dependence in LO or NLO means the large (small) contribution of high 
order diagrams, but it is not the case. The strong or weak scale dependence in LO or NLO indicates only 
the same for the higher orders, which nevertheless can be numerically small or large at some particular fixed 
scale. 
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assumed to be independent on the virtualities. This approximation considerably simphfies 
the calculations. 



The conventional NLO parton model approach with coUinear approximation works quite 
reasonably for one-particle distributions and the total cross sections; at the same time there 
is the serious disagreement with the data on azimuthal correlations and on the transverse 
momentum distributions of the heavy-quark pair. The reason is quite clear. Consider the 
processes of the fusion of two partons, having transverse momenta qix and g2T, into QQ pair 
with transverse momenta pir and p2T- In LO parton model qiT + (l2T = P1T+P2T = 0, and the 
distribution on the transverse momenta of the quark pair coinsides with the distribution on 
the total transverse momentum of the initial partons. NLO gives a correction, numerically 
not very large because the shapes of pt distribitions in LO and NLO are rather close to each 



other, see p, In this case the qualitative features in NLO should be the same as LO 
ones, so the final transverse momentum distributions of the heavy-quark pair only slightly 
differ from (5-functions, being in contradiction with the existing data 0. 

It was shown 0, that the conventional NLO parton model can describe the experimental 
data on azimuthal correlations and the momentum distributions of the heavy-quark pair, 
assuming the comparatively large intrinsic transverse momentum of incoming partons {kx 
kick) about 1 GeV/c. This procedure is, to a large extent, arbitrary. Let priQQ) be the 
total transverse momentum of the pair. For each event the heavy quark pair, QQ, is boosted 
firstly to the rest in the longitudinal direction. Then a second transverse boost is performed, 
giving the pair a transverse momentum priQQ) = ^i,T + fc2,T; ^i.r and A;2,r are the transverse 
momenta of the incoming partons, which are chosen randomly, with their moduli distributed 
according to 

Alternatively, one can proceed as in the previous case, but giving the additional transverse 
momentum ki^T + ^2,t to the whole final-state system (at the NLO in QCD, this means the 
QQ pair plus a light parton), and not to the QQ pair only. These two methods lead to 
similar results 0. 

However the fc^ kick method has no good theoretical background. The general phe- 
nomenological expression for the heavy quark production can be written^ as a convolution 
of the initial transverse momenta distributions, /(q'it) and /(g2T), with squared modulo of 
the matrix element, 

(^QCd{QQ) OC j d^qiTd^q2Tl{<ilT)I{<l2T)\M{qiT,q2T,PlT,P2T)? ■ (5) 



Here there are two possibilities: 

i) the typical gluon transverse momenta are much smaller than the transverse momenta of 
produced heavy quarks, qiT <C PiT, or 

ii) all transverse momenta are of the same order, qix ~ PiT- 

In the first case one can replace both initial distributions I {qiT) by 5-functions. It reduces 
^We omit here for the simphcity aU not important factors. 
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the expression (^) to the simple form (collinear approximation): 

(TcoldQQ) (X |M(0,0,piT,P2T)r (6) 

in total agreement with Weizsaecker- Williams approximation. 

In the second case we can not a priory expect good results from the Weizsaecker- Williams 
approximation. However it gives the reasonable numerical results in some cases. 

The kick ^ discussed above effectively accounts for the transverse momenta of incident 
partons. It is based on the expression which symbollically reads 

(TMckiQQ) oc /(giT)/(g2T) ® |M(0,0,pit,P2t)|' . (7) 

The main difference from the general QCD expression, Eq. (§), is that due to the absence of 
QiT in the matrix element the value (k^) in Eq. has to be differently taken for different 
processes (say, for heavy flavour production with comparatively small or large pt)- The 
reason is that the functions I{qiT) in general QCD expression decrease for large Qjfj^ as a 
weak power (see next Sect.), i.e. comparatively slowly, and the q^q. dependence of the matrix 
element is more important. 



3 Heavy quark production in 
kx-factorization approach 

Consider another approach, when the transverse momenta of the incident gluons in the 
small-x region result from the diffusion in the parton (gluon) evolution. The diffusion is 
described by the function (^(x, q^) giving the gluon distribution at the fixed fraction of the 
longitudinal momentum of the initial hadron, x, and the gluon virtuality, q^. It is approx- 
imately determined by the derivative of the usual structure function. Although generally 
it is the function of three variables, x, and g^, the transverse momentum dependence is 
comparatively weak since q^ ~ —q^ for small x in LLA in agreement with g^-dependences 
of structure function. Note that due to QCD scaling violation the value ip{x, q^) for the 
realistic structure functions increases with decreasing of x, and qx becomes important in the 
numerical calculations. 

The exact expression for g^ gluon distribution can be obtained as a solution of the 
evolution equation which, contrary to the parton model case, is nonlinear due to interactions 
between the partons in small x region. The calculations [2l| of the g^ gluon distribution in 



leading order using the BFKL theory result in difference from our <^(x, g^) function only 
about 10-15%. 

Here we deal with Eq. (|^) with the matrix element accounting for the gluon virtualities 
and polarizations. Since it is much more complicated than in the parton model we consider 
only LO of the subprocess gg QQ which gives the main contribution to the heavy quark 
production cross section at small x, see the diagrams in Fig. 2. The lower and upper ladder 
blocks present the two-dimensional gluon functions ip{xi^q\) and V9(x2,g|). In terms of the 
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conventional DGLAP parton distributions they can be determined as 

.2^ _ ^^^39[xG{x,q^)] 



(8) 



Such a definition of '^[x^ q^) makes it possible to treat correctly the effects arising from the 
gluon virtualities. 

Strictly speaking we have to include double logarithmic Sudakov-like form factor under 
the derivative in Eq. (8) [|20[. It was shown in that the probability to find a parton with 
a given longitudinal momentum fraction x and transverse momentum kt within the LLA 
may be written as 

d 



fa{x,kl) - 



(9) 



where a = g,q; T is the Double Logarithmic (DL) Sudakov form factor and 



^a[X,q 



4^2 7rVa(x,g') . 



(10) 



The first factor in (9) is evident. It corresponds to the real contribution/emission in the 
DGLAP evolution. Changing the scale from /x^ to //^ + one produces a new parton with 
k^ between /i^ and /i^ + So 

d 

fa{x,k1) = -^[xa{x,k1)] . (11) 

Next factor 

(C = Cp = {N^ — '^)/'2Nc for the quarks and C = Ca = Nc for the gluons) accounts for the 
virtual loop DGLAP contribution, which is needed to provide an appropriate normalization 
of the parton wave function and to satisfy the sum rules (i.e. the conservation of momentum, 
fiavour, etc.). T plays the role of survival probability, i.e. the probability not to emit an 
extra partons (gluons) with the transverse momenta q[ C [kt, fL\. 

Note that in DGLAP equation, written for the integrated (including all kt < fi) partons, 
there was a cancellation between the real and virtual soft gluon DL contributions. Emitting 
a soft gluon with the momentum fraction {1 — z) ^0 one does not change the x-distribution 
of parent partons. Thus the virtual and real contributions originated by the 1/(1 — z) 
singularity of a splitting function P{z) cancels each other. 

Contrary, in unintegrated case the emission of soft gluon (with > kt) alters the trans- 
verse momentum of parent (t-channel) parton. Eq. (9) includes for this effect through the 
derivative dT{kl, (j?) / dk'^. 

Thanks to this last derivative one obtains a positive unintegrated distribution fa{x,k^) 
even at rather large x, where (due to the virtual term in DGLAP) the integrated parton 
density a^x^jj?) decreases with /x^. 
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Unfortunately, the factor T is known with the DLog accuracy only. When kt > jj, there 
are no double logarithms at all. So for fcf > /i we put T = 1. 

In terms of Feynman diagrams in a physical (axial) gauge the factor T comes from the 
self-energy graphs. In Feynman gauge the self-energy contribution has no double logarithms. 
Such a double logarithms come from the graphs with the gluons which embrace the hard 
blob (an example is shown in Fig. lb). Any diagram with a soft gluon emitted from one 
external line of hard blob (say, gluon ^2) and absorbed on another external line has the same 
double logarithm, except of the colour factor. The sum of colour factors, corresponding to 
all three graphs (soft gluon is absorbed by the heavy quark, antiquark or gluon qi), is equal 
to the colour factor for self-energy diagram. This way one obtains in Feynman gauge the 
same expression (9). 

Thus the true probability to find a parton with the momentum fraction x and transverse 
momentum = k"^ is given by the Eq. (9). (see for more detail). To simplify the 



discussion we postpone the more exact Eq. (9) to the end of Section 4 and start with the 
Eq. (8) omitting the double logarithmic factor T. 

A different way to introduce the gluon qx distributions is suggested in Ref. ^6] and based 
on the Fourier transform of the structure functions. 

In what follows we use Sudakov decomposition for quark momenta pi^2 through the 
momenta of colliding hadrons pa and pb {pa = Pb — ^) transverse ones pi^2T- 

Pl,2 = Xi^2PB + yi,2PA + Pl,2T- (13) 

The differential cross section of heavy quarks hadroproduction has the form:0 

11/" 

j d'^^iTd'^^^rK^iT + q2T - Pit - P2t) 



dcTpp 



dyldyld?piT(Pp2T (27r) 

X ^^^{qly)^{qlx)\MQQ\\ (14) 

Here s = 2paPb , i?i,2T are the gluon transverse momenta and yl2 are the quark rapidities 
in the hadron-hadron c.m.s. frame, 

xi=^e-2^i, X2= ^e-y'^, x = Xi + X2 

y,= evt, 1/2 = ^ y = yi + y2 (15) 

2 2 2 

'^l,2T ~ "^Q + Pl,2T- 

|MqqP is the square of the matrix element for the heavy quark pair hadroproduction. 
In LLA kinematic 

qi ^ ypA + qiT, q2 ^ xpB + q2T- (16) 

so 

ql - -qlr^ ql - -qlr- (i7) 



^We put the argument of as to be equal to gluon virtuality, which is very close to the BLM scheme p7[|; 
(see also pi). 
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(The more accurate relations are qf = —j^, ~ working in the kinematics 

where x,y ^ 0). 

The matrix element M is calculated in the Born order of QCD without standart simpli- 
fications of the parton model. In the axial gauge p^A^ = the gluon propagator takes the 
form D^i,{q) = df,y{q)/q^, 

d,M = - (q'PB + q''p'B)/iPBq) . (18) 



For the gluons in t— channel the main contribution comes from the so called "nonsense" 
polarization ^f"^, which can be picked out by decomposing the numerator into longitudinal 
and transverse parts: 

= ^{p'bPa + P>b)/s + ^ 2p%p^Js ^ g;^ . (19) 

The other contributions are suppressed by the powers of s. It is easy to check that in axial 
gauge (^8|) Psd^u = and p^d^iu — —Qt/v- Thus we effectively get 



d,M ^ -2^. (20) 

ys 

Another way to obtain the same result is to use the transversality condition. Since the sum 
of the diagram Fig. 2a-c is gauge invariant the product qi^M^ = 0, here denotes the 
amplitude of the gluon qi and hadron interaction described by the lower part of Fig. 2a-c. 
Using the expansion ([16|) for the qi momentum we obtain 



which leads to Eq. (EDf) or 



y 



d,M ^2^, (21) 

xys 



if we do such a trick for the vector ps too Q Both these equations for d^j^u can be used but 
for the form (^) one has to modify the gluon vertex slightly (to account for several ways of 
gluon emission - see ref. [Q): 



xys 



(xys - ql^) qiT " qItQ2T + 2a; {qiTq2T) Pb ■ (22) 



As a result the colliding gluons can be treated as aligned ones and their polarization vectors 
are directed along the transverse momenta. Ultimately, the nontrivial azimuthal correlations 
must arise between the transverse momenta piT and p2T of the heavy quarks. 

From the formal point of view there is a danger to loose the gauge invariance in dealing 
with the off mass shell gluons. Say, in the covariant Feynman gauge the new graphs (similar 

^Having in mind this trick one can say that the matrix element is gauge-invariant in fcy-factorization 
approach. However the polarization vectors of incoming gluons qi and q2 are not arbitrary taken but fixed 
as —Qit/v and — gjr/^^i respectively (see for more detail). 
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to the " bremsstruhlung" from the initial or final quark line, as it is shown in Fig. 2d) may con- 
tribute in the central plato rapidity region. However this is not the fact. Within the "semi- 
hard" accuracy, when the function <^(x, q^) collects the terms of the form aj(lng^)"(ln(l/x))™' 
with n + m > k, the triple gluon vertex (^) includes effectively all the leading logarithmic 
contributions of the Fig. 2d type [25, 18|. For instance, the upper part of the graph shown in 



Fig. 2d corresponds in terms of the BFKL equation to the t-channel gluon reggeization. Thus 
the final expression is gauge invariant (except a small, non- logarithmic, 0{as) corrections) 

5 

Although the situation considered here seems to be quite opposite to the parton model 



there is a certain limit |^, in which our formulae can be transformed into parton model 
ones, namely when the quark transverse momenta, pi,2T; are much larger than the gluon 
ones, gi or- 



4 Results of numerical calculations 

Eq. (|H|) enables to calculate straightforwardly all distributions concerning one-particle or 
pair production. To illustrate the difference between our approach and the conventional 
parton model we present first of all the results of calculations of charm production {rric = 



1.4 GeV |2|, 0) with high pr at three energies, y/s = 1 TeV, 10 TeV and 10^ TeV and with 



the same value of 

= ^ = 0.02 , (23) 
i.e. pt = 10 GeV, 100 GeV and 10'^ Gev for the above energies. 



We will use QCD scales fip = ttit = ynT-l + Pt ^^id fipt = "^c, i-e. fixed coupling with Nf 
= 3 and A = 248 MeV. 

However there exists a problem coming from the infrared region, because the functions 
</9(x,g|) and <^{y,qi) are unknown at the small values of g^2- Moreover, for the realistic 
gluon structure functions the value (f{x, q^) is positive in the small-x region and negative in 
the large-x region. The boundary between two regions, where ^p{x^ q^) =0, depends on the 
g^, therefore the relative contributions of these regions is determined by the characteristic 
scale of the reaction, say, by the transverse momentum pt- 

To avoid this additional problem, we present the numerical calculations both in the kx- 
factorization approach and in the LO parton model, using the "toy" gluon distribution 

xG{x,q^) = {l-xfHq^/Ql) , (24) 

for g2 > Ql and xG{x, q^) = for q^ < Ql with Ql = 1 GeV\ 

After the calculations of charm production cross sections using Eqs. (1) and (8), (14) one 
can see, that the values da/dxT at xt = 0.02 are about 4-5 times larger in the fc^-factorization 
approach than in the LO parton model. Really this difference should not be considered as 



^Taking the gluon polarization vector as — ^^y/y we can completely eliminate the leading logarithm terms 
arising from Fig. 2d. 
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very large because, as it was discussed, an essential part of NLO and NNLO corrections 
is already included in the /c^-factorization, and it is known that the sum of LO and NLO 



contributions in the parton model is about 2-3 times larger than the LO contribution 



To show the most important values of the variables gi,2T in Eq. {\L4^ , as well as the 
kinematical region producing the main difference with the conventional parton model, we 
plot by dots in Fig. 3 the results of the fcr-factorization approach with the restrictions 
\qi,2t\ < Qmax, as a function of Qmax- The running coupling is fixed as as{ml) instead of 
C(s{qi2) in Eq. The dashed lines in Fig. 3 are the conventional parton model predictions 



with QCD scales fip = yf^c +Pt-> ^^d = rric- One can see that the /cT-factorization 
predictions exceed the parton model results when qmax > Pt- 

Let us check now that the /cj'-factorization predictions really coincide with the parton 
model ones, when the quark momenta, Pi^2T, are much larger than the gluon ones, gi,2T 



28[| . However we have to compare them at the same values of structure functions. When 
we submit Eq. (0) to the conditions |gi,2r| < Qmax and neglect the gi^2T dependence in the 
matrix element, we recover the parton model approximation, and get the result proportional 



to xG{x, ql^ax) ' y^hji Qmax) dircct consequence of Eq. (|) 

xGix, q^) = xG{x, Ql) + (^(x, ql)dql . (25) 

4v27r'^ JQl 



At the same time the original parton model yields the result proportional to xG{x,fi'jp) ■ 
yG{y, fi'jp) with fip = \Jm1 + p^. Hence we expect the parton model to coincide with our 
calculations for the gluon distribution Eq. (|2^) , ^ "^c and |gi,2T| < qmax after multiplying 
by an appropriate factor: 



da 
dxr 



da 

dxr 



PM 




(26) 



The right-hand side of Eq. presented in Fig. 3 by solid curves shows a good agreement 
with /cy-factorization approach (open dots) even when qmax is only slightly smaller (at the 
highest energy) than px- The same values da/dxp, as in Fig. 3, but differentiated with 
respect to In q^ax are presented in Fig. 4 for /cy-factorization approach. It exhibits, especially 
at the high energies, the logarithmic growth with qmax, until the value qmax ~ Pt- There is 
the narrow peak in this region, which provides about 70-80 % of the da/dxx cross section 
integrated over the whole gi,2T phase space. Its physical nature seems to be quite clear. 
There are two kinematical regions for t-channel and w-channel diagrams. Fig. 2 a,b, giving 
comparatively large contributions to da/dxT for the high energies and high pt heavy quark 
production. One of them comes from the conventional parton model kinematics when both 
initial gluon transverse momenta, gi,2T; are very small compared to pi,2T; see Fig. 5a. Another 
large contribution appears from the region where, say, qiT ~ Pit, whereas q2T and p2T 
are comparatively small, as it is shown in Fig. 5b. In this case the quark propagator, 
l/(pi — gi) —mq) = l/{q2 —P2) —^q), is close to the mass shell and gives rise to the narrow 
peak shown in Fig. 4. The general smallness of such processes comes from high-virtuality 
gluon propagator in Fig. 5b, and it is of the same order as in the case of Fig. 5a, even with 
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Table 1: The ratios of cc pair production in /cr-factorization approach and in LO parton 
model 



yi, TeV 


0.3 


1 


10 


100 


1000 


Rtot 


4.0 


4.0 


4.0 


3.9 


3.9 


R{xt = 0.02) 


3.4 


4.5 


5.5 


5.4 


5.2 



some suppression of the diagram Fig. 5a, where the contribution from the region of large 
rapidity difference of two produced heavy quark is suppressed by the fermion propagator. 

The diagram shown in Fig. 5b can be considered |]l| as one of the most important contri- 
bution to the NLO parton model in the high energy limit, because it have spin-one gluon in 
the i-channel. Due to all these factors, combinatorics and the interference between diagrams, 
the summary contribution of the peak to the total da / dxx value is several times larger than 
the LO parton model contribution. In the Table 1 we present the calculated ratios of the 
total heavy quark pair production cross section, Rtot, and the ratio of da/dxr, R{xt) at xt 
= 0.02 (integrated over rapidity). 

One can see that the relative contribution of the discussed peaks firstly increases due 
to increase the phase space. At the energy high enough {^/s ~ 10 TeV this contribution is 
saturated. After that the relative contribution of the LO parton model increases logarith- 
mically with and have to dominate at the extremely high energies and pxi in academical 
asymptotic. 

These results are confirmed by those presented in Fig. 6, where we reproduce the data 
from Fig. 3 for LO parton model and fcy-factorization approach with the condition |(3'i,2t| < 
Qmax- For comparison we show by stars the fcr-factorization predictions for all values qiT 
with the only condition |g2T| < <imaxi versus qmax- The last results are above the LO parton 
model even at not too large g2r because the peaks, discussed above, are included into the 
integral over qi^. 

Let us note that the calculation of da/ dxx at a/s = 10 TeV, xt = 0.02, and with restriction 
|?i,2r| > Pt/2 (see Fig. 5c) gives only about 2% of the total value of da/dxx- 

It is necessary to note that the essential values of qi^2T increase in our calculations when 
the transverse momentum, pt of the detected c-quark increases. At initial energy high 
enough, gi,2T ~ Pt- In the kr kick language it means that the value also increases. 

The values da/dpx are presented in Fig. 7. Here the lower dashed curves correspond to 
the LO parton model, solid curves — to the fc-r-factorization approach with unintegrated 
gluon function (p{x, q^) given by the simplified expression (8) while the dottted and dash- 
dotted curves embody the double-logarithmic T-form factor (12) and use Eqs. (9) and (10) 
(see the discussion in section 3). Recall that from the diagram viewpoint the T-factor sums 
up the diagrams containing in the Feynman gauge the gluon lines embracing the hard blob as 
it is shown in Fig. lb. In the axial gauge the corresponding Double Log contribution comes 
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from the parton (gluon) self-energy insertion, i.e. from the term proportional to 6{1 — z) 
in the DGLAP evolution equation. The T-factor reproduces a large part of the virtual 
NLO corrections with respect to the conventional PM and diminishes the cross section. It 
suppresses the contribution coming from small values of parton virtuahties, ^12^ "^r = 
rriQ + in Eq. ([T^) . Such a suppression becomes more important at very large pt as in 
this case the PM contribution, which is collected from the logarithmically large interval of 

9i,2 ^ Pt (/q2 d.11,2/11,2)} starts to dominate. 

Recall that before we have no " scale" problem within the fcj-factorization approach since 
we deal with the explicit integral (|1^ over the parton virtualities g^2- The dependence of 
QCD matrix element M{qiT,q2T,PiT,P2T) on the gi^2T values provides the convergency to 
the integral "automatically" fixing the "scale" ^p. 

After T-factor is "switched on" one gets the problem of scale back. The scale fi^ is 
not fixed in ( |l2l) at LO accuracy. The reasonable /i^ values are expected to be somewhere 
between m'^ and s = M^-^ = xys. The numerical difference between these values is very 
small, as one can see from Fig. 7. 

The rapidity distributions of prodiced charm quark is presented in Fig. 8. Here again 
the dashed curves are the results of LO parton model calculations, the sohd curves show the 
results of the fcr-factorization approach with unintegrated gluon function ip{x, q^) given by 
the simplified expression (8), while the dotted curves embody the double- logarithmic T-form 
factor (|12D with T(m|,) and use Eqs. (9) and (10). One can see that the last form factor 
decrease the cross section calculated in the /cT-factorization approach about 1.5 times. 

5 Conclusion 

We have compared the conventional LO Parton Model (PM) and the /cT-factorization ap- 
proach for heavy quark hadropro duct ion. In order to concentrate on the specific features of 
the approach and not to deal with the pecularities of the parton distributions a simple gluon 
density ( P^ has been chosen. 

It has been shown that the contribution from the domain with strong ordering 
{.<ii,2T -C rriT = \l"i^Q + Pt) coincides in /cT-factorization approach with the LO PM predic- 
tion. However a very numerically large contribution comes besides this in /cT^-factorization 
from the region q\^2T > "^t- It is kinematically related to the events where the transverse 
momentum of heavy quark Q is balanced not by the momentum of antiquark Q but by the 
momentum of the nearest gluon. 

This configuration is associated with NLO (or even NNLO, if both q\;xT > ^t) corrections 
in terms of the PM with fixed number of flavours, i.e. without the heavy quarks in the 
evolution. Indeed, as it was mentioned in |1]], up to 80% of the whole NLO cross section 
originates from the events where the heavy quark transverse momentum is balanced by the 
nearest gluon jet. Thus the large "NLO" contribution, especially at large pt, is explained 
by the fact that the virtuality of the t-channel (or w-channel) quark becomes small in the 
region q^- ~ px and the singularity of the quark propagator l/(p — ^) — ^^q) in the "hard" 
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QCD matrix element, M{qiT,q2T,PiT,P2T), reveals itself. 

Including the double logarithmic Sudakov-type form factor T into the definition of un- 
integrated parton density one takes into account an important part of the virtual loop 
NLO (with respect to PM) corrections. Thus we demonstrate that /cT-factorization collects 
already at LO the major part of the contributions which play the role of the NLO (and even 
NNLO) corrections to the conventional PM. Therefore we hope that the higher order (in as) 
correction to the fcr-factorization would be rather small. 

Another advantage of this approach is that a non-zero transverse momentum of QQ- 
system {pit + P2T = ?ir + i?2r) is naturally achieved in fcr-factorization. The typical value of 
this momentum (fcr-kick) depends on the parton structure functions/densities. It increases 
with the initial energy (/c^-kick increases with the decreasing of the momentum fractions x, y 
carried by the incoming partons) and with the transverse momenta of heavy quarks, pt- Thus 
one gets a possibility to describe a non-trivial azimuthal correlation without introducing a 
large "phenomenological" intrinsic transverse momentum of the partons. 

A more detail study of the heavy quark production (including the correlations) in the kr- 
factorization approach based on the modern realistic parton distributions will be published 
elsewhere. 
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Figure captions 



Fig. 1. Heavy quark production in the LO parton model (a) and in fcr-factorization 
approacli with Sudakov T-factor (see discussions in Sect. 3). 

Fig. 2. Low order QCD diagrams for heavy quark production in pp [pp) collisions via 
gluon-gluon fusion (a-c) and the diagram (d) formally violating the factorization, that is 
restored within logarithmic accuracy. 

Fig. 3. The values of da/dx^ for charm hadropro duct ion at xt = 0.02 in /cT^-factorization 
approach as a function of upper limits of integrals in Eq. (p^ ) (points); the conventional 
parton model values (dashed lines) and the values of right-hand side of Eq. ( p6D (solid 
curves) . 

Fig. 4. The values of da /dxx/d In qmax for charm hadroproduction at xt = 0.02 in 
/c^-factorization approach as a function of upper limits of integrals in Eq. (|I^ ) (points). 

Fig. 5. The diagrams which are important in the case of one-particle distributions of 
heavy quark with large pt- The situation similar to the LO parton model (a). The case, 
possible in NLO parton model, large pt of the quark is compensated by hard gluon and 
the fermion propagator is near to mass shell (b). The numerically small contribution, when 
large transverse momenta of heavy quarks are compensated by two hard gluons, whereas the 
fermion propagator is near to mass shell (c). 

Fig. 6. The same as in Fig. 3 (points and dashed curves) together with the values of 
da/dxT (stars) for the case when only g2T upper limit integration is bounded by qmax- 

Fig. 7. pT-distributions of c-quarks produced at different energies. Dashed curves are the 
results of LO parton model. Solid curves are calculated with unintegrated gluon distribution 
(p{x,q^) given by Eq. (10), dash-dotted and dotted curves from Eqs. (11) and (12) for ij? 
values in Eq. (12) equal to s/4 and m^, respectively. 

Fig. 8. Rapidity distributions of c-quarks produced at different energies. Dashed curves 
are the results of LO parton model. Solid curves are calculated with unintegrated gluon 
distribution (p{x,q'^) given by Eq. (10) and dotted curves from Eqs. (11) and (12) for /i^ 
values in Eq. (12) equal to m^. 
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